Discovery of exceptional points in the Bose-Einstein condensation of gases with 

attractive 1/r-interaction 
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The extended Gross-Pitaevskii equation for the Bose-Einstein condensation of gases with attrac- 
tive 1/r-interaction has a second solution which is born together with the ground state in a tangent 
bifurcation. At the bifurcation point both states coalesce, i.e., the energies and the wave functions 
are identical. We investigate the bifurcation point in the context of exceptional points, a phe- 
nomenon known for linear non-Hermitian Hamiltonians. We point out that the mean field energy, 
the chemical potential, and the wave functions show the same behavior as an exceptional point in a 
linear, non-symmetric system. The analysis of the analytically continued Gross-Pitaevskii equation 
reveals complex waves at negative scattering lengths below the tangent bifurcation. These solutions 
are interpreted as a decay of the condensate caused by an absorbing potential. 

PACS numbers: 03.75.Hh, 02.30.-f, 34.20.Cf, 04.40.-b 



I. INTRODUCTION 

It is a well known property of Bose-Einstein conden- 
sates with attractive interatomic interactions that sta- 
tionary solutions to the Gross-Pitaevskii equation exist 
only in certain regions of the parameter space governing 
the physics of the condensates. For example, for the case 
of an attractive s-wave contact interaction it was theo- 
retically predicted (see, e.g., [l|, H, Q), and experimen- 
tally confirmed [1, S Sli ^^at the condensate collapses 
when, for given negative scattering length, the number 
of particles becomes too large. In an alternative exper- 
iment p] the collapse was induced by tuning the scat- 
tering length in the vicinity of Feshbach resonances b 
adjusting an external magnetic field. Huepe et al. [T 
have shown that the critical parameter values where col- 
lapse occurs in fact correspond to bifurcation points of 
the solutions to the stationary Gross-Pitaevskii equation 
and analyzed the linear stability of the solutions. Quite 
recently it was demonstrated [1(| that similar behavior 
persists in Bose-Einstein condensates where, in addition 
to the short-range (van der Waals-like) interaction a long- 
range "gravity-like", attractive 1/r interaction is present. 
Such monopolar quantum gases could be realized accord- 
ing to O'Dell et al. [ll[ by a combination of 6 appro- 
priately arranged "triads" of intense off-resonant laser 
beams. There, too, when crossing the borderlines in the 
parameter space, spanned by particle number, scattering 
length and trap frequency, at which a stationary ground 
state of the extended Gross-Pitaevskii equation comes 
into being, a second, excited, state of the condensate is 
born in a tangent bifurcation. 

It is the purpose of this paper to examine these bifur- 
cations from the point of view of the theory of "excep- 
tional points" , since at the bifurcation points both the 
eigenvalues and the wave functions of the two states are 
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identical, a situation well known from studies of excep- 
tional points 121, which have been investigated theoret- 
icallyElH lljll [13, [3 [S ll and experimentally 
pll [22 . 23, 24, 2511 in a wide variety of physical systems. 

Exceptional points can appear in systems described 
by non-Hermitian matrices which depend on a two- 
dimensional parameter space. At critical points in the 
parameter space (exceptional points) a coalescence of two 
eigenstates can occur, where both the eigenvalues and the 
eigenvectors of the two states pass through a branch point 
singularity and become identical. There is only one lin- 
early independent eigenvector of the two states at the 
exceptional point. 

In the quantum mechanics of linear Schrodinger equa- 
tions, exceptional points are investigated in open systems 
in which resonances exist. Examples are discussed, e.g ., 
for complex atoms in laser fields 26] , a double S well [23] , 
the scatte ring of a beam of particles by a double barrier 
potential [2g], and the hydrogen atom in static external 
fields [2^ . The existence of resonances is important for 
the occurrence of exceptional points in quantum systems 
because the coalescence of two eigenstates is not possible 
in the case of Hermitian Hamiltonians describing bound 
states. In the latter case, there is always a set of or- 
thogonal eigenstates which never become identical. The 
situation is different for resonances, which can, e.g., be 
described by now-Hermitian Hamiltonians. In this case, 
the eigenstates have not to be orthogonal and a coales- 
cence can occur. 

In this paper we will reveal the existence of excep- 
tional points also in quantum systems described by non- 
linear Schrodinger equations. As a model system we 
choose Bose-Einstein condensates with attractive 1/r- 
interaction and concentrate on the case of self trapping, 
i.e., condensation without external trap, which is a fea- 
ture of such systems [l3, [HI • It is shown that the bifur- 
cations of the two stationary solutions to the nonlinear 
Gross-Pitaevskii equation at critical physical parameter 
values exhibit the typical structure known from stud- 
ies of exceptional points in linear systems. We should 
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point out that the effect is not restricted to condensates 
with 1 / r- interaction. Our model system has the great ad- 
vantage that simple approximate analytic solutions exist, 
which can be used for the investigation of the exceptional 
points. The branch point singularity structure can be 
seen directly from the analytic terms we obtain, which is 
only possible for self-trapped condensates without exter- 
nal harmonic trap. 

The Gross-Pitaevskii equation for self-trapped Bose- 
Einstcin condensates with monopolar l/r- interaction is 
introduced in Sec. [Ill In Sec. Illli we give an overview on 
exceptional points in linear systems and their properties, 
which can be used to identify them. The complex vicin- 
ity of the bifurcation point is investigated in Sec. IIVI to 
reveal the branch point singularity structure of the bifur- 
cation and to identify it as a "nonlinear version" of an 
exceptional point. The physical interpretation of a com- 
plex absorbing potential in the Gross-Pitaevskii equation 
at scattering lengths below the critical value at the tan- 
gent bifurcation is discussed in Sec. |Vl Conclusions are 
drawn in Sec. IVII 



II. 



BOSE-EINSTEIN CONDENSATES WITH 
1 /r-INTERACTION 



In this section we briefly review the equations and re- 
sults for self-trapped Bose-Einstein condensates with at- 
tractive 1 / r- interaction which are necessary for the sub- 
sequent analysis of exceptional points. The extended 
Gross-Pitaevskii equation without external trap poten- 
tial reads 



\ CLu J |r - r' 




where the natural "atomic" units introduced in [10| were 
used. Lengths are measured in units of a "Bohr radius" 
a„ and energies in units of a "Rydberg energy" which 
are given by 



E,, — 



2ma,, 



respectively, where u determines the strength of the 
atom-atom-interaction 11], and m the mass of one bo- 
son. In Eq. ([1]), e is the chemical potential, a the s-wave 
scattering length and N the number of bosons. As was 
pointed out in jlQ] with the use of the scaling property 
of the system, the physics of a self-trapped condensate 
is only determined by one parameter, namely iV^a/a„, 
which will be necessary for the identification of the ex- 
ceptional point. A further quantity necessary for our 
discussions is the mean field energy of the self-trapped 



condensate, which reads 



E[iIj\ ^ N I d^r 7/'*(r) ( -A^ + inN— \tP{r) 



"/^■vM(r:i^),<r,. (2) 



The solutions we are interested in, viz. the two states 
emerging at the tangent bifurcation, are radially symmet- 
ric. One of them is the ground state. Thus we concen- 
trate on radially symmetric solutions of Eq. Q . In this 
case both analytic calculations via a variational method 
and numerically exact computations can be carried out. 



A. Variational solutions 

An analytic approximation for the two wave functions 
of the condensate which emerge at the tangent bifurca- 
tion can be obtained using a variational principle. Cal- 
culations with a Gaussian type orbital were performed 
by O'Dell et al. |11| an d compared with numerically ac- 
curate solutions in [lO|. The trial wave function is given 
by 



ijj{r) 



where 
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3/2 



(3) 
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is the normalization constant and the variation is per- 
formed with respect to the width represented by k. The 
two stationary values. 
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of the mean field energy ^ are obtained at 
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ill Gu 



(6) 



represents the variational approximation for the 
ground state of the condensate, the second solution (ex- 
cited state) is labeled E-. Furthermore, the analytical 
expressions for the chemical potentials of the two solu- 
tions are given by 



£±_ ^ _± 



5±4jH- :f 



(7) 



respectively. The tangent bifurcation is obvious from the 
equations ^ and ([7]). For negative scattering lengths 
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FIG. 1: (Color online) Chemical potentials of the two so- 
lutions. The results of the variational and the numeri- 
cal exact calculations are shown. Two solutions emerge in 
a tangent bifurcation at a critical value of the parameter 
N'^a/uu- In the variational approximation the critical value 
is N^a/uu = —3tv/8, whereas the numerical exact solutions 
emerge at N^a/au = -1.0251 The "atomic" units intro- 
duced in 1101 are used. 



to fmax ~ 15 is sufBcient. The wave function obtained 
with the initial condition ^'(fo) = 1 is not normahzed, 
which is indicated with the tilde. Because of the scaling 
property [ic| of the problem, the normalization of a nu- 
merically obtained solution ^ is possible with the help of 
a normalization factor v. As can be seen with a simple 
calculation, the transformation 



(12) 



leaves the system ([5]) of differential equations invari- 
ant. With the normalization condition ~ N and 
the scaling invariance (|12p . the normalization factor is 
given by 



= ll^-ll^ = 47r / |*(f)pf2df 



(13) 



and the properly scaled and normalized wave function 
^'(r) is obtained as 



with the scaled scattering length 



(14) 



with N^a/ttu < — 37r/8 there is no (real) result for E and 
e, at the critical value N'^a/a^ = — Stt/S both solutions 
have the same value (£'+ = e+ = e_), and above 
— 37r/8, we obtain two different real solutions. The vari- 
ational chemical potential is shown as a function of the 
scattering length parameter N'^a/a^ in Fig. [TJ 



B. Exact calculations 

For the search of radially symmetric wave functions, 
the radial part of the integro-differential equation ^ is 
written in the form of two differential equations. 



*"(r) + -*'(r) = -U{r)-^{r) + 87r6|^' (r)p*(r) , (8) 
r 



U"(r) + -U'(r) = -87r|^'(r)|^ 
r 

with the parameter 

an 



(9) 



(10) 



The A^-dependence in ([T]) has been absorbed in the wave 
function = ^/N->p. The two second order differential 
equations can be transformed to first order equations and 
integrated, e.g., with a Runge-Kutta-Merson algorithm 
starting at, e.g., vq = 10^^^ slightly greater than zero 
with the initial conditions 

*(fo) = 1 , *'(ro) - , U{fo) = uo , U'{fo) = . (11) 

The value uq must be determined such that '^(f) van- 
ishes in the limit f oo. Numerically, integration up 



(15) 



An explicit calculation of the integral (|13p is not required. 
Using the asymptotic behavior of the numerically com- 
puted potential [/(f) for large radial coordinates 



~ 1 

£+2—- , 
1^ r 



the factor u/N can be calculated from 

f .. -2 

N 



lim 



(16) 



(17) 



The correctly scaled chemical potential can also be de- 
termined with the help of the approximation (|16p : 



e 

iV2 



(^)'.l™ (;7(f)+fJ7'(f)) . (18) 



Only for the mean field energy a further integral is re- 
quired. Using the virial theorem [l], [ll[ we can obtain it 
as 



E 



(£1 

Af2 



= -47r 



■f{r)fU{f)r''df. (19) 



The numerical exact chemical potential is also shown in 
Fig. [T] It agrees, qualitatively, with the result of the 
variational approach, however, the tangent bifurcation is 
shifted to a higher critical scattering length N'^a/au — 
-1.0251 

We will show in Sec. lIVI that in both the variational ap- 
proach and the exact calculations the tangent bifurcation 
at the critical scattering length is a branch point singu- 
larity of the wave functions, i.e., an exceptional point. 
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III. EXCEPTIONAL POINTS 

To make our presentation self-contained we recapit- 
ulate the essential concepts on exceptional points and 
their properties known from linear systems, which will 
be important for the comparison with the coalescence 
of two quantum states in a nonlinear Schrodinger equa- 
tion in Sec. IIVI Usually complex symmetric matrices or 
complex symmetric matrix representations of Hamiltoni- 
ans are used to investi gate exceptional points in quan- 
tum systems (see, e.g.. [l5l I27I [291]'). Here, we treat the 
more general case of non-symmetric matrices, considered 
before [s^, HH, H^l also for systems with time reversal 
symmetry breaking [H, [s^, HH] , and compare the phase 
behavior of the eigenvectors with the symmetric case. It 
is instructive to discuss a simple two-dimensional model 
which exhibits exceptional points. The example is pro- 
vided by the 2 x 2-matrix 



(a) 



M(A) = 



1 A 

+ A -1 



(20) 



In (j20p . A is a complex parameter, and the complex num- 
ber c is introduced to distinguish between the general 
non-symmetric case (c ^ 0) and a symmetric matrix 
(c = 0) in the discussion below. The eigenvalues read 



3i,2(A) =±Vl + cA + A2 



(21) 



and are obviously two branches of one analytic function. 
Corresponding non-normalized eigenvectors of the two 
eigenvalues are given by 

-A 



a;i.2(A) 



1 T VI + cA + A2 



(22) 



There are two exceptional points for Xa,b = — c/2 ± 
\/ (? 1^ — 1. At these parameter values, both the eigenval- 
ues 61^2 and the eigenvectors X\-i pass through a branch 
point singularity. For the further discussions we concen- 
trate on the case \ — \a- Then the degenerate eigenval- 
ues have the value ei,2(AA) =0 and the two eigenvectors 
X\ 2 are 



a;i,2(Ay 



c/2-vA^ 
1 



(23) 



An important property of exceptional points, which fol- 
lows from the branch point singularity structure, is the 
permutation of the two eigenvalues if the exceptional 
point is encircled in the parameter space [l2|. To illus- 
trate this effect we follow the paths of the eigenvalues 61^2 
in the complex plane for a circle with radius q around the 
critical value Aa, which can be represented by 



Ag((^) = -c/2+ VcV4-l + ee''^ 



(24) 



An approximation for small radii q ^ \2^J / A — 1 1 leads 
to 



3i(A,(^)) = yW^^7I^e'(^/2+V4) ^ 
32(A,(^)) = yW?7^^e'(^/2+5V4) _ 



(25) 
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FIG. 2: (Color online) Eigenvalues of the two-dimensional 
model defined in Eq. (|2Up with c = 1 for a circle around the 
exceptional point, (a) Circle in the parameter space with 
radius q = 10"'^. The starting point A(0) is marked with an 
open square and the direction of progression is indicated with 
the arrow, (b) Path of the eigenvalues 61,2 for the parameter 
values on the circle. The eigenvalues which belong to the first 
parameter value A(0) are labeled with open squares and the 
arrows point in the direction of progression. 



If a full circle in the complex parameter space A is tra- 
versed, the paths of the eigenvalues form a semicircle as 
can be seen from Eq. (pS)) . Neither of the two eigenval- 
ues passes through a closed loop. The first arrives, after 
the parameter space loop, at the starting point of the 
second one and vice versa. Two circles in the A-space 
{lp — Q . . . 47r) are required to obtain a full circle of one 
of the eigenvalues. The situation is shown in Fig. O for 
c—\. In the parameter space the circle (l24l) is traversed 
once with a radius g = 10~^, as plotted in Fig. [2] (a). 
The paths of the two eigenvalues in the complex e-space 
is shown in Fig. [5] (b). The semicircle structure expected 
from approximation (|25p is clearly visible. 

For exceptional points in linear systems there is a 
noteworthy difference for complex symmetric and non- 
symmetric matrices. In the first case, which is obtained 
for c = in the model (j20p and which was realized in the 
resonances investigated in microwave cavities [22| and 
in the hydrogen atom in crossed electric and magnetic 
fields [231 , there is always a distinct phase behavior of the 
eigenvectors. If an exceptional point is encircled, the two 
eigenvectors are interchanged (as expected from the be- 
havior of the energy eigenvalues) and, additionally, one 
of the two eigenvectors changes its sign. This effect is 
demonstrated in [Tsj and can be summarized with, e.g.. 



circle 



(26) 



Four circles around the exceptional point in the param- 
eter space are needed to restore the original situation 
[xi,X2\. A direct verification of the effect was performed 
with microwave cavity experiments (2^ , in which a visu- 
alization of the wave functions was possible. The phase 
behavior can be obtained for the vectors Xi normalized 
with the Euclidean norm without complex conjugation 
Ni = yj xfxi, where xf is the transpose of the vector Xi 
(cf. |20(|). This normalization fixes the phase and is the 
method required for the comparison with the nonlinear 
system studied in this paper (cf. Sec. IIVI and appendix). 
For the non-symmetric case the change in sign does not 
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FIG. 3: (Color online) Phase of the product pi2 for a circle in 
the parameter space of type (1241) with radius g = 10~^ around 
the exceptional point, (a) Symmetric matrix with c — 0. The 
phase changes by n indicating the change in sign of one of 
the eigenvectors from (|26p . (b) Non-symmetric matrix with 
c = 1. There is no change in sign as noted in (|27|l . 



appear for a vector whose phase has been fixed with the 
same method, i.e., the Euchdean norm without complex 
conjugation. A non-symmetric matrix is obtained in the 
model (j20p for c ^ 0. Then, there is only a permutation 
of the eigenvectors, which can be summarized with 



[Xl,X2\ 



[X2,xi] 



(27) 



This is shown analytically for a small circle in the ap- 
pendix. Note that the phase behavior (p7|) does not rep- 
resent the "geometric phase" of an eigenvector during 
the circle around an exceptional point, which has been 
calculated in Ref. [3l| . 

A demonstration of the phase behavior is possible with 
the product 



Pl2 



1 

1 



X2 



(28) 



The non-diagonal matrix ensures, that the product pi2 
does not vanish in the symmetric case. If the change in 
sign is present, the phase of pi2 changes its value by tt 
when a full circle around the exceptional point is tra- 
versed. In Fig. [3] the behavior of the product is shown 
for the example (j20p . The change in sign is obvious in 
the symmetric case c = and does not appear for the 
non-symmetric choice c — 1. 



IV. ANALYTIC CONTINUATION OF THE 
GROSS-PITAEVSKII SYSTEM 

In the extended Gross-Pitaevskii equation ([T]), both 
the variational and the numerical calculations suggest the 
existence of an exceptional point at the tangent bifurca- 
tion. For the critical parameter value, we obtain two 



identical states. The chemical potentials, the mean field 
energies and the wave functions are identical. In con- 
trast to exceptional points in linear systems described in 
Sec. mil the coalescence of the two states of the nonlin- 
ear Schrodinger equation ([1]) can appear for a purely real 
wave function at a real energy in a one-dimensional pa- 
rameter space. This is a consequence of the nonlinearity 
of the Schrodinger equation. 

If we want to check whether or not the degeneracy in 
the nonlinear system has the same branch point singu- 
larity structure as exceptional points in open quantum 
systems described by a linear Schrodinger equation, we 
have to extend the parameter N'^a/uu to complex values 
and to investigate the complex vicinity of the degeneracy. 

The analytic continuation of the Gross-Pitaevskii sys- 
tem is a nontrivial task. The Gross-Pitaevskii equation 
([!]) contains the square modulus of the wave function tjj 
and is therefore a non- analytic function of ip. It has been 
argued that the tempting simple replacement of 1-01^ with 
■0^ is valid only with the assumption that the entire wave 
function is real valued [s^, [33|. Here, we suggest the fol- 
lowing procedure for complex wave functions. 

Any complex wave function can be written as 



V'(r) 



(29) 



where the real functions a{r) and (3{r) determine the 
amplitude and phase of the wave function, respectively. 
The complex conjugate and the square modulus of ijj{r) 
read 



V'*(r) = e"M-''3W , |V'(r)|^ 



(30) 



Using the ansatz ([29)1 the Gross-Pitaevskii system can be 
written as two coupled nonlinear differential equation for 
a{r) and (3{r), however, without any complex conjugate 
or square modulus. These equations can now be contin- 
ued analytically by allowing for complex valued functions 
a{r) and f3{r). This implies that Eq. ([50]) without com- 
plex conjugate of a(r) and /3{r) is formally used for the 
calculation of tjj* and |'(/'(r)p, and thus the square modu- 
lus of ip can become complex. The physical interpretation 
of a complex absorbing potential in the Gross-Pitaevskii 
system will be discussed below in Sec. |Vl 

It should be noted that the above procedure is analo- 
gous to the complex scaling method, viz. the replacement 
r — > re'^ in the Hamiltonian, used for the calculation of 
resonances in open quantum systems [ssl . [39j . For states 
on the left hand side of operators there is no complex con- 
jugation of the re'^ arising from the scale transformation 
but complex conjugation is applied to the intrinsically 
complex part of the function. 

In the following we will encircle the branch point sin- 
gularity at the tangent bifurcation. For the variational 
approach discussed in Sec. Ill Al the Eqs. (O, ([6]), and 
([7|) are straightforwardly extended to complex values of 
the parameter N'^a/am yielding complex results for the 
chemical potential, the mean field energy, and the wave 
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functions. The situation is more complicated for the ex- 
act calculations in Sec. lIIBI because the differential equa- 
tions ([8]) and ([9]) do not directly depend on the scaled 
scattering length A^^a/a„ but on the parameter b. For 
complex b the differential equations ([8]) and ([9]) can be 
continued analytically resulting in complex wave func- 
tions. Then, the two complex parameters uq and b are 
determined in a multidimensional root search problem, 
to fulfill the condition of vanishing wave function in the 
limit r oo, and to achieve the given value of A^^a/a„. 



A. Branch point singularity structure of the 
energies 

The variational solutions ([5]) for the mean field energy, 
and ([7]) for the chemical potential show a branch point 
singularity at the bifurcation point. For a circle around 
the critical value of the complex parameter N'^a/au the 
typical permutation of the two solutions is expected and 
found, as is shown in Fig. |4] for the variational approx- 
imation as well as for the numerical exact calculations. 
The circle in the parameter space is defined by 



(^ = 0...27r, (31) 



where (A^^a/o„)c is set to the critical value 
(7V^a/atj)var = — 37r/8 = —1.1780 in the variational 
and to {N'^a/au)num ~ —1.0251 in the numerical exact 
calculations, respectively. We used the radius g — 10^^. 
Both cases are shown in Fig. |4] (a) and (b). In Fig. 
S] (c) and (d) the variational and numerical solutions 
for the chemical potential are shown. The results were 
obtained for different parameter values located on the 
circle. As can be seen in the figure, the permutation of 
the two values of the chemical potential appears very 
clearly. Each of the two solutions e+ and e_ traverses a 
path in the complex energy plane similar to a semicircle. 
A fractional power series expansion of the variational 
result (O for small radii. 
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FIG. 4: (Color online) Chemical potential and mean field 
energy for a circle in the complex N'^a/uu parameter space. 
The first row shows the parameter space circle of type (|3ip 
with g = used for the variational (a) and the numerical 
exact (b) calculations. The starting point is marked with an 
open square and the direction of progression is indicated with 
an arrow. In the second row, the chemical potential e/N'^ for 
the variational (c) and numerical (d) solutions is plotted. The 
variational (e) and numerical (f) mean field energies E/N^ 
are presented in the third row. Each of the two solutions 
is drawn with a different line. The results obtained for the 
first parameter value on the circle are marked with an open 
square and the arrows point in the direction of progression. 
A permutation of the two solutions is present for the mean 
field energy as well as for the chemical potential indicating 
the existence of an exceptional point. All energies are given 
in units of the "Rydberg energy" Eu mentioned in Sec. IIIl 
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confirms this finding. The term • y^e'*^/^, which dom- 
inates the path of the eigenvalue for a value g ^ 1, leads 
to a semicircle. 

For the mean field energy the permutation is also 
present but the structure of the paths of the two solu- 
tions is different. As can be seen from the fractional 
power series expansion of the analytic solution ([5|) , 
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4 32 
97r ~ 97r2 



.y^%W2)i^ + 0(y^'*) , (33) 



the first order term with the phase term e^'^^'^ vanishes. 
The lowest non- vanishing order has the phase factor e"^, 
which leads to a closed curve for a complete circle in the 
parameter space (ip — 0...27r). The third order term 
is the lowest order responsible for a permutation. As 
a consequence, it becomes more and more difficult to 
see a permutation for decreasing radii g. Nevertheless, 
the permutation of the two values is present. The same 
result can be seen in Fig. |4] (e) and (f). The dominating 
structure is given by the circle following from the second 
order term but the permutation is clearly visible even for 
the small radius g = 10^^. 

The exceptional point discussed here was discovered for 
a real value of the parameter N'^a/uu, however, it is still 
an isolated point in the two-dimensional parameter space 
spanned by the real and imaginary parts of N'^a/au- This 
finding is in full agreement with non-Hermitian linear 
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FIG. 5: (Color online) Phase arg(Oi2) of the overlap integral 
Oi2 defined in Eq. (|34|l . The path in the parameter space 
N'^a/au is a circle with the critical value as center point and 
Q — 10~^. The angle on the circle is denoted by if. Both the 
variational and the numerical results show that the phase of 
Oi2 returns to its initial value at the end of the circle and 
that the change in sign from (|26p is not present as expected 
(see text) . 



systems, in which the co-dimension of exceptional points 
is two (see, e.g., [lOl). The analytic continuation is re- 
quired to confirm the branch point singularity structure. 
As mentioned above, the appearance on the real axis is 
possible because of the nonlinearity of the Schrodinger 
equation. 

B. Behavior of the wave functions 

The eigenf unctions of the extended stationary Gross- 
Pitaevskii equation ([T]) are not orthogonal. They behave 
more like the eigenstates of a non-symmetric linear sys- 
tem and the phase behavior (j27p without a change in sign 
is expected. A possibility to check this is to calculate the 
complex overlap integral for the two non-orthogonal nor- 
malized states 

Oi2 = 47r / Mr)Mr)r''dr (34) 
Jo 

for parameter values located on a circle of the type ([M)) . 
For this calculation it is important to note that the nu- 
merical wave functions calculated with the method de- 
scribed above have a phase, which is fixed by the nor- 
malization factor V determined by Eq. (|13p without com- 
plex conjugation for the analytic continuation of the ex- 
tended Gross-Pitaevskii model (cf. Eq. ((30)l ). A similar 
procedure is used for the variational result, where the 
normalization constant A in is calculated with a nor- 
malization integral without complex conjugation. 

If one of the two states changes its sign during the 
traversal of the loop, the phase of the complex value O12 
changes its value by tt similarly to the phase behavior of 
the product pi2 from Eq. (pS)) . If this is not the case, the 
phase returns to its original value at the end of the loop. 

In Fig. [5] the phase of the integral ([M)) is plotted for 
the circle defined in Eq. ((3T|) with radius g — 10~^. The 
result shows clearly that the phase of O12 returns to its 
initial value after one circle around the exceptional point, 
demonstrating that no change in sign of the eigenvectors 



occurs. Thus, the wave functions behave in the same way 
as for exceptional points in linear systems described by 
a non-symmetric matrix (see Sec. IIII[) . 

V. DECAY OF THE CONDENSATE 

The Gross-Pitaevskii equation ([1]) depends on the scat- 
tering length N'^a/au of the contact potential. In Sec. 
llVl that parameter has been extended to complex values, 
although the physical scattering length is always real. 
There is, however, a physical situation where complex 
continuation is needed for real TV^a/a.^, viz. when the 
scattering length is below the critical value at the tangent 
bifurcation. For N^a/uu < — 37r/8 = —1.1780 the varia- 
tional solutions ^ and ([7]) for the mean field energy and 
the chemical potential become complex. Those parame- 
ters become complex also for the analytically continued 
exact calculations at iV^a/a.^ < —1.0251. A complex 
chemical potential as a signature of a decaying conden- 
sate has already been discussed in Ref. [40| . 

The standard physical interpretation of complex 
eigenenergies is that they describe decaying resonances 
in open systems or systems with absorbing potentials. 
The imaginary part of the energy E is related to the 
width r, decay rate A, and lifetime T of the resonance by 
T — hX — h/T — — 2Im E. Both the real and complex 
wave function of the Bose-Einstein condensate above and 
below the tangent bifurcation vanishes at large radius r, 
i.e., the wave function cannot describe a decay of the con- 
densate by outgoing particles. How is a decay possible? 
The solution is that for a complex wave function ^(r), 
using the rules ([50)1 for complex conjugate and square 
modulus, the effective potential 

Kff (r) = 8^7V-|V'(r)p -2N [ d'r' ^fi^ (35) 
flu J \r~r'\ 

becomes an absorbing potential. Thus, the decay is an in- 
ternal collapse of the condensate. The physical interpre- 
tation might be that at large negative scattering lengths 
the contact potential is so attractive that the atoms of 
the condensate form molecules or clusters as discussed, 
e.g., in Refs. 0,110) El- The decreasing number of atoms 
means the decay of the condensate. The imaginary part 
of the chemical potential as a function of the scattering 
length is presented in Fig. [S] for the variational as well as 
for the exact calculation. As was the case for parameter 
values above the critical point, there are two solutions. 
Below the critical N'^a/uu one is the complex conjugate 
of the other. 



VI. CONCLUSION 

We have investigated the complex vicinity of the 
tangent bifurcation appearing in the extended Gross- 
Pitaevskii equation for self-trapped Bose-Einstein con- 
densates with attractive 1/r interaction. The chemical 
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FIG. 6: (Color online) Imaginary part of the chemical po- 
tential for real N^a/uu below the bifurcation point obtained 
with the variational and the numerical calculations. Both the 
continuation of the ground state and the continuation of the 
excited state are plotted. One is the complex conjugate of the 
other. 



potential and the mean field energy show the typical be- 
havior of energies at an exceptional point. Both quanti- 
ties pass through a branch point singularity at the criti- 
cal value of the system's parameter 7V^a/a„. An analysis 
of the phase behavior of the corresponding wave func- 
tions showed that the bifurcation point has the same 
properties as an exceptional point of a non-Hermitian 
non-symmetric matrix. Thus, we conclude that we have 
identified the bifurcation point as a "nonlinear version" 
of an exceptional point. 

For the investigation of the bifurcation point an an- 
alytic continuation of the Gross-Pitaevskii system was 
used. As we have discussed, an analytic continuation of 
the wave functions is not only required for complex val- 
ues of N'^a/uu but is necessary for real parameter values 
below the bifurcation point. Then, a complex absorbing 
effective potential is obtained indicating an internal col- 
lapse of the potential. An investigation of the dynamics 
of the system similar to Refs. [1, [9| will give more insight 
in the stability of the states and is the topic of current 
studies. 

In this paper we have identified the bifurcation points 
appearing in self-trapped monopolar Bose-Einstein con- 
densates as exceptional points. There is good reason to 
believe that quite generally the critical parameter values 
of attractive Bose-Einstein condensates where collapse of 
the condensates sets in are associated with exceptional 
points. The existence of a bifurcation point is known for 
condensates in a harmonic trap without l/r- interaction 
[1, [1] and for condensates with dipole-dipole interaction 
and harmonic trap. We have checked that the bifurcation 
points appearing in variational approximations of these 
systems are exceptional points. Clearly this should be 
confirmed by detailed analysis including numerically ex- 
act calculations for attractive Bose-Einstein condensates 
based on the techniques applied in this paper. 

Acknowledgments 



graduiertenforderung of the Land Baden- Wiirttemberg. 



APPENDIX: PHASE OF THE EIGENVECTORS 
ON THE CIRCLE 

To be able to compare the phase behavior of the wave 
functions of the nonlinear Gross-Pitaevskii system with 
the eigenvectors of a matrix, we fix the phase of the vec- 
tors with the same method as is discussed in Sec. llVlfor 
the wave functions of the analytic continuation of the ex- 
tended Gross-Pitaevskii model. This is done with t he Eu- 
clidean norm without complex conjugation Ni — \J xfxi 
of the eigenvectors, which is the adequate counterpart of 
the normalization factors v from Eq. (|13p and A from Eq. 
(|31), where (formally) the complex conjugation is ignored. 
Then, the normalized eigenvectors have the form 



a;i,2(A) 



1 



A^i,: 



-A 



1 T Vl + cA + A2 



with 



A^i.: 



+ cA + A2 



For a small circle around the exceptional point A^ = 
— c/2-f ^c2/4 — 1 we calculate the fractional power series 
expansion in p^/^ of the normalized vectors and obtain 
in the non-symmetric case (c ^ 0) 



K/rj^e'^fi/^ + 0(£i3/4) 
K?7V^e'^/2 ^ Q(^3/4) 

^yge''^/^ ^ 0(^3/4) 




with K — ■y/c2/4 — 1, and rj ~ — 2ck. 

The traversal of a single circle {ip = . . . 2n) obviously 
leads to a permutation of the eigenvectors as summarized 
in (I26l) . 

In the symmetric case (c = 0), a fractional power se- 
ries expansion of the two normalized eigenvectors looks 
different. The result is 



.19^/8 






+ 


25/4 






ei5-/8 


23/" 


el/4 + 


25/4 


gillT/8 


e-i'.'/o 


ei5x/8 


23/4 


el/4 + 


25/4 




e-i*./4 


e-/s 


23/0 


el/4 -t 


25/4 



pl/4eiW4 + 0(f?3/4)^ 
^gl/4eW4 + 0(^3/4)^ 



With a circle around the exceptional point ((p = . . . 2ti) 
it can be seen directly that Xi{2'it) — — a;2(0), and 
X2{2tt) — a;i(0), i.e., the permutation rule ([^5]) for a 
symmetric matrix is reproduced in this picture. 
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